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ABSTRACT 

Context. The aromatic infrared bands (AIBs) observed in the mid infrared spectrum of galactic and extragalactic sources are attributed 
to Polycyclic Aromatic Hydrocarbons (PAHs). Recently, two new approaches have been developed to analyze the variations of AIBs 
in terms of chemical evolution of PAH species: Blind Signal Separation (BSS) and the NASA Ames PAH IR Spectroscopic Database 
fitting tool. 

Aims. We aim to study AIBs in a Photo-Dissociation Region (PDR) since in these regions, as the radiation environment changes, the 
evolution of AIBs are observed. 

Methods. We observe the NGC 7023-North West (NW) PDR in the mid-infrared (10 - 19.5 fim) using the Infrared Spectrometer 
(IRS), on board Spitzer, in the high-resolution, short wavelength mode. Clear variations are observed in the spectra, most notably the 
ratio of the 1 1 .0 to 11.2 fjm bands, the peak position of the 11.2 and 12.0 /jm bands, and the degree of asymmetry of the 11.2 fim band. 
The observed variations appear to change as a function of position within the PDR. We aim to explain these variations by a change 
in the abundances of the emitting components of the PDR. A Blind Signal Separation (BSS) method, i.e. a Non-Negative Matrix 
Factorization algorithm is applied to separate the observed spectrum into components. Using the NASA Ames PAH IR Spectroscopic 
Database, these extracted signals are fit. The observed signals alone were also fit using the database and these components are com- 
pared to the BSS components. 

Results. Three component signals were extracted from the observation using BSS. We attribute the three signals to ionized PAHs, 
neutral PAHs, and Very Small Grains (VSGs). The fit of the BSS extracted spectra with the PAH database further confirms the attri- 
bution to PAH^ and PAH" and provides confidence in both methods for producing reliable results. 

Conclusions. The 1 1.0 yum feature is attributed to PAH+ while the 1 1.2 /im band is attributed to PAH". The VSG signal shows a char- 
acteristically asymmetric broad feature at 1 1.3 yum with an extended red wing. By combining the NASA Ames PAH IR Spectroscopic 
Database fit with the BSS method, the independent results of each method can be confirmed and some limitations of each method are 
overcome. 
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1. Introduction 

Polycyclic Aromatic Hydrocarbons (PAHs) are carbonaceous 
macromolecules which were postulated to be present in the in- 
terstellar medium (ISM) in the 1980s (jLeger & P uget 1984[ 
AUamandola et al.|[T985j [Puget & Leger||1989[ [AllamandoTal 
et al. II 1989 1 and have since undergone an intense investigation 
in astronomy. The state of the art and recent activity in the field 
of interstellar PAHs is well illustrated by the book "PAHs and 
the Universe", IJoblin & Tielens] (^ll ). Astronomical PAHs are 
generally considered to contain roughly 50 - 1 00 C atoms and 
have an abundance of a few 10"^ per H atom( Tie lens 2008 l 



Because of their nanometer size, the absorption of one far- 
ultraviolet (FUV) photon is sufficient to heat PAH molecules 
to high temperatures causing them to emit characteristic bands 
called Aromatic Infrared Bands (AIBs) which peak near 3.3, 
6.2, 7.7, 8.6, and 1 1.2 //m ( Tielens|2008 and references therein). 



* This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. 



PAHs are abundantly present in the diffuse ISM, reflection nebu- 
lae (RNe), planetary nebulae, protoplanetary disks, and galaxies. 
Observations of PAHs in Photo-Dissociation Regions (PDRs, 
[Hollenbach & Tielens 199 9), which are transition regions be- 
tween atomic and molecular gas, but still strongly affected by the 
FUV photons, are of particular interest. The UV flux decreases 
when moving from the neutral atomic gas to the dense molecu- 
lar cloud and PAH populations will also evolve as the UV flux 
changes ( jSloan et al.|1999; Joblin et al.|199 6). This effect is best 
studied in the mid-infrared (5-15 fim), where PAHs emit most 
strongly. 



Each AIB is charachteristic of a PAH vibrational mode (e.g. 
Puget & Leger[T989 , AUamandola et al.|1985[|Hony et al.|2001|l. 



the 3.3, 8.6, and 10 - 15 jum features are due to the C-H stretch- 
ing, in-plane and out-of-plane bending modes while the 6.2 and 
7.7 jum features are mainly due to the C-C stretching modes. 
These features have been observed to show strong vaiiation in 
peak position, width of band (FWHM), and symmetry (Peeters 
et al. 20021. The changing ratio of the 8.6 and 11.3 fim fea- 
tures was discovered first and attributed to a change in the rela- 
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live abundances of neutral and ionized PAHs ( |JobUn et al.|1996[ 



Sloan et al. 1999 1. Soon after, the 7.7/11.3 jim ratio was observed 



to vary as well, which was also attributed to the charge state 
of the PAHs. Later, using Infrared Space Observatory (ISO), 
Peeters et al. ( 2002[ l catalogued the variations of the main fea- 
tures in the 6-9 /im range and empirically divided them into 
groups based on specific spectral properties. It was found that 
each of these groups was representative of certain classes of ob- 
jects: Class A included HII regions, RNe, galaxies, and non- 
isolated Herbig stars. Class B included isolated Herbig stars, 
PNe, and two post AGB stars, while class C included only two 
post AGB stars. Observations of the 10-15 fim region have been 
analyzed by Hony et al. (2001 ) in terms of the solo, duo, trio, 
or quartet out-of-plane (OOP) bending mode of either PAH", 
PAH^, or some combination of the two. While the position and 
profile of these bands are quite characteristic the relative inten- 
sities do vary a lot, indicating variations in the edge structure of 
the aromatic molecules. 



Analysis and interpretation of astronomical observations is 
supported by dedicated laboratory studies and quantum chemi- 
cal studies. These studies are being carried out in many groups 
around the world, each using a different technique (see Oomens 
( 201 l| l for a review of laboratory and experimental studies). In 
most cases, the absorption is studied at a low temperature in 
an inert matrix. There are also some gas-phase experiments that 
have been carried out at higher temperatures. These experimen- 
tal studies have been extended by quantum chemical calculations 
using Density Functional Theory (DFT), to species not accessi- 
ble in laboratory studies. The DFT approach is used to determine 
the frequencies and intensities of vibrational modes. Recently, 
these models have been used to calculate spectra for PAHs from 
54 to 130 C atoms ( IBauschlicheTeFal. 2008 2009). This size 
range is particularly relevant for comparison to observations 
of space-based PAHs. An extensive database has been created 
by the Astrophysics and Astrochemistry Laboratory at NASA 
Ames, which includes mid-IR to far-IR spectra of many different 
PAHs including large molecules, varied levels of ionization, and 
irregular shapes ( Bauschlicher et al.|2010[ l. Malloci et al. ( 2007| l 
used time dependent DFT methods and quantum-chemical cal- 
culations to report computed molecular properties of PAH emis- 
sion for 40 molecules, available on an online database. Mulas| 
|et al.| ( p006a) then modeled the PAH emission, which give po- 
sitions and intensities of specific PAHs in different radiation en- 
vironments. The band profiles of some PAH emission were also 
calculated by Mulas et al. (2006b). 



Recently, the improved sensitivity of Spitzer Space 
Telescope has brought a wealth of observations of AIB features. 
NGC 7023 is a well studied and bright IR source where PAH 
variations are kn own to occur ([Cesarsky et al. 1996| l. Werner 
et al. (2004 1 and Sellgren et al.| (2007i used Spitzer' s Infrared 
Spectrograph (IRS) in the Short-High (SH), Short-Low (SL), and 
High-Low (HL) modes to further observe the full mid-infrared 
range of NGC 7023. They observed all the classical AIBs above 
5 /im in addition to finding new, weak emission features at 6.7, 
10.1, 15.8, 17.4, and 19.0 fim. ^erae et al.| ( |200T :) observed NGC 
7023, along with three other PDRs, using Spitzer's IRS-SL. The 
spectra were analyzed using a class of methods called Blind 
Signal Separation (BSS), which identifies elementary spectra 
from spectral cubes. Using BSS on spectra of NGC 7023-NW, 
three component signals were recognized, PAH cations, neutral 
PAHs, and a third carrier which Berne et al.l ilOQl ) attributed 



what the exact nature of VSGs are, it has been proposed that 
they could be PAH clusters (Rap acioli et al.|2()05 i. 

This paper presents a study of the PDR NGC 7023-NW, 
which aims to put observational constraints on the origins of the 
profiles and variations of the 10 - 15 jim spectra. This study will 
complement and restrict previous results from quantum chemi- 
cal calculations, ISO spectroscopy observations, Spitzer IRS ob- 
servations, and PAH models. We analyze high resolution data 
from Spitzer's IRS-SH ( [Werner et al.|[2004| , at a resolution, 
R - A/AA - 600, using the NASA Ames PAH IR Spectroscopic 
Database and a BSS method to separate the emitting components 
of the PDR. 

After briefly outlining the observational methods in Section 

2, the Blind Signal Separation method is described in Section 

3, including the application of BSS to the data (Section 3.2). 
Section 4 presents our main results and compares these with pre- 
vious studies of the region. Section 5 provides a comparison to 
the fit with the NASA Ames PAH IR Spectroscopic Database. 
Next, in Section 6, we discuss the implications of our results and 
propose strong candidates to explain the spectral variations of 
the 10-15 /im region. Section 7 discusses briefly the nature of 
the VSG carrier and section 8 gives our concluding remarks. 



2. Observations 

NGC 7023-NW is a PDR located 40" to the northwest of the ex- 
citing star, HD 200775, seen in Figure[T] HD 200775 is a magni- 
tude 7 Herbig Be star and is located 430 pc from the sun ( van den 
|Anckeretal.|I997| l. There are 3 PDRs in NGC 7023 located east, 
south and northwest of the exciting star The northwest PDR is 
the brightest of the 3 PDRs. 

NGC-7023-NW was chosen for this study in view of the in- 
teresting results of the analysis of low resolution data of this re- 
gion performed by Berne et al. ( 2007| l who have shown that the 
AIB spectrum can be separated into 3 main components. To fur- 
ther this study, high-resolution short wavelength (IRS=SH, 10 - 
19.5 yum) observations were obtained from the Spitzer Heritage 
Archive (SHA). The data was reduced using the CUbe Builder 
for IRS Spectra Maps (CUBISM), provi ded by NASA's Spitzer 
Space Telescope tools ( [Smith et al.|2007| l. 



to evaporating Very Small Grains (VSG), following earlier as- 
signment by [Rapacioli et al!] ( [2005| l. Although it is yet unclear 



3. Methods 

3.1. Blind Signal Separation 

Blind Signal Separation (BSS) is a class of methods used in 
many scientific fields to separate source signals from observed 
linear combinations of these signals e.g. separating brainwaves 
or unmixing recordings in acoustics. This has only recently been 
applied to astronomy (Nuzillard & Bijaoui 2000 Berne et al. 
2007). In the case of observing a PDR, the resultant spectrum 
at any spatial point is a superposition of all elements emitting 
in the designated wavelength range and BSS can separate these 
components. There are three main methods to perform BSS: 
Independent Component Analysis (ICA), Non-negative Matrix 
Factorization (NMF), and Sparse Component Analysis (SCA). 
Based on the results of Berne et al. (2007) and the added spectral 
resolution of the data, the NMF method was selected to perform 
the analysis. If X is a matrix containing the observed spectra, as- 
suming that each spectrum is the result of a linear combination 
of source signals, we can write X ^ WH, where H is the matrix 
of source signals, and W is the matrix of mixing coefficients. 
The goal of NMF is to recover W and H based on X only. The 
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Fig, 1, IRAC 8jum image of NGC 7023 ( [Werner et al.||2004| l. 
Highlighted with a white box is the IRS-SH field of view around 
the PDR NGC 7023-NW. The star, HD 200775, is marked by the 
green label. 



method we apply here is the same as Berne et al. (2007i where 
all details can be found. 



3.2. Application to NGC 7023 

Our astronomical data is comprised of 750 spectra, each taken 
at a different spatial location, and each spectrum including 869 
wavelength points. Figure |2] shows two spectra, before they are 
decomposed with BSS. There is a clearly defined separation be- 
tween the 11.0 and 11.2 yum emission features as well as distinct 
features at 12.0, 12.7, 13.5, and 14.1 fim with additional fea- 
tures at longer wavelengths. Among the features at longer wave- 
lengths are the PAH 16.4 yum feature, the blended PAH and Ceo 
feature at 17.5 fim, and the pure Ceo feature at 18.9 jum, which 
are further discussed in Sellgren et al. (2010 1. The relative inten- 
sities of the H2 lines vary with position in the nebula, represent- 
ing a non-linear component in our spectra. Another non-linear 
aspect of the spectra is the onset of the dust grain continuum 
caused by heating from the source star. These non-linear compo- 
nents cannot be analyzed by BSS methods, since the method de- 
mands a linear combination of signals. Therefore, we have cho- 
sen to exclude the emission at wavelengths greater than 15 /im 
where the continuum from dust grains is present. We have also 
clipped the H2 lines everywhere in the cube by hand and replaced 
them by a linear interpolation. The resulting spectra were then 
analyzed using the NMF algorithm from |Lee & Seung ( 2001| l 
with both divergence and Euclidian distance optimization, see 
|Bern6 et alT| ( |2007i ) for details. 

We investigated the possibility of 3, 4, 5, and 6 component 
source signals. Irrespective of the particular minimization tech- 
nique, when attempting to separate 4, 5 and 6 sources, there are 
2 or more signals which are very similar (linear combinations of 
each other) and at least one signal that is pure noise. Therefore, 
we can conclude that there are 3 significantly different spectral 
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Fig. 2. Observed spectra in pixel position [15,28] (blue) and 
[23,17] (red), containing H2 lines and thermal continuum. 




Fig. 4. Spectra taken at three random spatial positions (left). The 
solid line represents the original observed spectra while the red 
overlapping circles (thick red line) represent a linear combina- 
tion of the BSS extracted spectra. On the right is the matching 
residuals for each plot. 



components responsible for the AIBs in NGC 7023-NW. This 
result confirms the findings of Berne et al. ( 20Q7\ that there are 
only 3 source signals in this PDR. 



4. Results of Blind Signal Separation 

4.1. Extracted Source Signals 

The final extracted spectra are shown in Figure |3] To increase 
confidence in these results, and ensure that this solution is not 
a random local minimum, the same analysis was repeated 100 
times using different random initializations. These 100 spectra 
shared the same general shape, but varied in intensity, especially 
in the 11.0 - 11.3 fim region. The average spectra of the 100 it- 
erations is plotted with a red line in Figure [3] and will be used 
for the remainder of our analysis as the final BSS extracted sig- 
nals (H matrix). We can also estimate the error at each point in 
the spectrum using these results (Figure|3]l. The BSS method has 
the most difficulty separating the signals in the 11.0 to 11.3 //m 
range due to the strong changing spectral gradients there. This 
results in large errors in this range (see Appendix |A] for discus- 
sion on unmixing artifacts). Since X » WH, we can estimate W 
by minimizing \\X - WH\\ using a standard least squares mini- 
mization. Figure |4] compares the observations in X and the final 
reconstruction of these observations with WxH. The reconstruc- 
tion is in good agreement with the observations. 

Using the weighting factors that come as a resultant matrix 
of the above reconstruciton (W), we can map the spatial distri- 
bution of each source signal separately (Figure |3]l. The spatial 
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Fig. 3. Bottom Panels: Extracted spectra using NMF, normalized at 11.2 jum. The vertical line represents the peak position of the 
estimated PAH° spectra. The red line represents the average spectra out of 100 iterations. The grey envelope shows the minimum 
and maximum spectra and the black envelope shows the l-cr error of the 100 iterations. Upper Panels: Spatial distributions of the 
weighting factors obtained by Least Squares Fitting of the observed spectra in the datacube using the BSS extracted spectra shown 
in red in the lower panel. 



distribution shows clear variation for the three emitting compo- 
nents. Signal 1 is most abundant in the middle of the PDR. Signal 
2 has its highest concentration closest to the source star (located 
at the bottom left of this image) and Signal 3 appears to trace the 
edge of the PDR farthest from the star The well defined regions 
where each signal is most concentrated implies a physical cause 
and gives further confidence that these results are not random. 



4.2. Carriers of the Extracted Spectra 

In this section, we will compare our results to the results of 
the low resolution study of the same region (Berne et al.'2007'; 
Berne and et al. 2010) to gain insight about the three extracted 



signals. Creating spatial contours of intensity for each signal al- 
lows us to compare the spatial distribution of our signals to the 
distribution of the three signals from the study of Berne et al.| 
(j2007) of the 5-15 fim low resolution spectra. The contours 
are created from the IRS-SH spatial distribution maps (Figure|3]l 
and overlaid with the spatial distributions (represented in color) 
of the IRS-SL results (Figure|5]). The three signals extracted here 
show a strong spatial correlation to the PAH^, PAH*\ and VSG 
maps o f|Berne et a l. (2007). The spatial distribution and the re- 
sults of Berne et al., (2007 ) seem to suggest that Signal 2 traces 
the distribution of PAH cations, Signal 1 the neutral PAH dis- 
tribution, and Signal 3 the distribution of VSGs. Although the 
spatial distributions of Signal 1, Signal 2, and Signal 3 corre - 
late well with PAH°, PAH+, and VSGs of [Berne et aL] ( |2007| , 
there are some small discrepancies, in particular, for the PAH" 



map. As discussed in Berne and et al. (2010 1, the degradation of 
spatial or spectral resolution always implies a loss in the quality 
of the NMF efficiency. Since Berne et al. (2007) have a higher 
spatial resolution, while here we have a higher spectral resolu- 
tion, none of the data-sets can be considered "better" and small 
discrepancies between the results of NMF are expected. 

Figure [6] co mpares the low-resolution source signal spectra 
of Berne et al. (2007i to the high-resolution source spectra ob- 



tained here (Figure[3]). The low resolution extracted spectra share 
all the major features with the high resolution spectra, specifi- 
cally the broad 11.3 fim emission feature in the VSGs, the 11.2 
IJ.m and 12.7 fim emission features in the neutral PAHs, and 
the 11.0 yum and broad 12.7 fim emission features in the PAH 
cations. Although evidence for the 11.0 fim feature was present 
in the IRS-SL observations of |Berne and et a L]( |2010[ ), it was 
not immediately attributed to PAH cations since the resolution 
was not high enough to fully resolve and separate the 11.0 and 
1 1.2 jum features. The IRS-SH spectra has not only validated the 
previous results of the IRS-SL observations but given us addi- 
tional spectral detail, to infer more about each signal and how it 
contributes to the observed spectra. This is of particular interest 
in the 11 fim region where the 11.0 fim and 11.2 fim bands are 
well resolved and isolated. We also observe a consistent broad- 
ening of the 11.3 feature in the VSG spectra. One main differ- 
ence in the spectra is the presence of an 11.0 fim satellite feature 
in Signal 1 that is not found in the IRS-SL PAH*' signal, which 
we suggest is an artifact caused by the inherent limitation of un- 
mixing and is compensated for by the sharp drop at 11.0 jum 
of Signal 3. The strong correlation between both the spectra and 
spatial distributions allow us to confidently identify Signal 1 as a 
pah" signal. Signal 2 as a PAH^, and Signal 3 as a genuine VSG 
spectrum. This also matches the physical description of the PDR 
with the ionized species closest to the source star 

5. Comparison with Database Fitting Analysis 

The BSS-method allows the separation of the observed spec- 
tra into three mathematically distinct components, without tak- 
ing the actual physical or spectroscopic properties of aromatic 
species into account. To further explore the carriers of the 
three signals resulting from the BSS analysis, we turned to the 
NASA A mes PAH IR Spectroscopic Database (BauschHcher] 
et al. 2010), which contains over 600 theoretical and experimen- 
tal PAH spectra. The existing extensive database of PAH spectra 
allows us to approach the analysis of the astronomical data from 
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(c) 



Beme et al. 



(2007b. (b) The contours 



Fig. 5. (a) The contours of Signal 2, overlaid with the distribution of PAH cations created froml] 

of Signal 1, overlaid with the distribution of neutral PAHs created from Berne et al. (2007 1, (c) The contours of Signal 3, overlaid 
with the distribution of VSGs created from| Berne et al.| ([2007). 
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Fig. 6. The top spectra are the results from|Berne and et al. (20101 using low-resolution IRS data, the bottom row of spectra are the 



current results, using high-resolution IRS. The vertical dashed black line indicates the 1 1.2 /im line position. 



a different perspective. Specifically, it allows us to link observa- 
tional properties of the infrared emission features to the molec- 
ular characteristics of the carriers. To that end, we first fit each 
BSS extracted spectrum with the PAH database. Second, we fo- 
cus on one observed spectrum from the IRS-SH data cube and 
use the database to fit this spectrum. In interpreting the results 
from this database analysis, it should be kept in mind that vi- 
brational modes in the mid-IR spectral range are characteristic 
for molecular groups and are not very sensitive to individual 
molecules. Hence, the goal of this database analysis is to pro- 
vide insights in trends rather than specific molecular identifica- 
tions. The trends of interest here involved the effects of charge, 
size, molecular geometry and symmetry - including the degree 
of compactness of the PAH families - on the different spectral 
components. As a corollary, the completeness of the database is 
therefore of lesser concern, as long as the relevant classes are 
well represented. 

The BSS method fits are "blind" in the sense that there is 
no a priori information about the nature of the signals built into 
the method. On the other hand, the fits of the database are based 
on spectra of actual aromatic molecules in specific charge states, 
structures, and sizes, allowing for a more direct interpretation of 
the results. 

The BSS results separated 3 mathematically distinct emit- 
ting signals and based on comparison with the signals extracted 
from Berne et aL] (2007 1, we have attributed the three signals to 
PAH^, PAH", and VSGs. However, there is no aspect of the BSS 
method that actually identifies the signals as ionized or neutral 
PAH species. By fitting these signals with the database, we can 



obtain an independent attribution of these distinct signals as PAH 
classes. 



5. 1 . Fit Parameters 

The database, at all versions, and the AmesPAHdblDLSuite can 
be obtained from www.astrochem.org/pahdb. Here, version 1.11 
of the theoretical component of the database and the November 
10, 2010 version of the IDL suite was used. Briefly, the spectra in 
the theoretical database correspond to about 600 different PAHs, 
ranging in size from C9H7 to C130H28. Note that the database is 
biased towards smaller species, with PAHs containing over 50 
carbon atoms making up roughly 24%. The database includes 
PAHs at different charge states (i.e. cations, anions and neutral 
species) as well as different symmetries of the same molecule. 
The fit included all C-H PAHs as well as polycyclic aromatic 
nitrogen heterocycles (PANH's; PAHs with one or more nitro- 
gen atoms substituted into their carbon skeleton), since Hudgins| 
et al.^ (2005j suggested that at least 1.2% of the cosmic nitro- 
gen is tied up in PAH molecules. The fit excludes PAHs with 
Oxygen, Magnesium, or Iron, where we note that no specific 
spectral evidence for the existence of such species has been 
found yet. 

Two approximations are made when fitting "observed" spec- 
tra with the database. First, the spectra in the database refer to 
absorption spectra at K, while the observed spectra are emitted 
by "hot" species. Due to anharmonicity, emission bands are ob- 
served to shift to the red with increasing temperature ( [Cherchneff 
|eraL][T992l |Cook & Saykally|[T998l ). Systematic experimental 
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Database 


BSS Components 


Percentage 


PAH+ 


PAH" 


VSG 


Large PAHs 


43 


48 


15 


Small PAHs 


57 


52 


85 


PAH Cations 


78 


15 


34 


PAH Neutrals 


16 


80 


21 


PAH Anions 


6 


5 


45 



Table 1. The numerical result of the fit using the PAH Database 
to fit the component spectra extracted using the BSS method. 
The percentage amounts of large (C > 50), small (C < 50), 
cation, neutral, and anion species are displayed for the PAH^, 
PAH*^, and VSG component signals. 



and quantum chemical studies on a very limited set of PAHs 
show that this redshift depends on the mode under considera- 



tion, the molecular structure and the temperature ( jJobUn et al. 
[1995) lOomens et aT]|2003) [Basire et al.|[20T0| . The out-of- 
plane bending modes are observed to shift by about 15 to 20 
cm ' for the small PAHs, pyrene and coronene between and 
900K and we will adopt this value for all PAHs in the database. 
Second, the database tools allow the user to specify a Gaussian 
or a Lorentzian profile with an assumed linewidth. However, the 
intrinsic line profile of PAH emission bands is distinctly non- 
Gaussian and non-Lorentzian due to the effects of anharmonic- 
ity and the accompanying red-shading ( jBarker et al. 1987; Pech 
|et al.||2002] l. Here, we will adopt Lorentzian profiles, fully re- 
alizing that this implies that this procedure will force the fit of 
the observed broader and red-shaded bands to blends by emis- 
sion from multiple species. We will adopt an intrinsic Lorentzian 
linewidth of 6 cm ' for the out-of -plane bending modes and note 
that this is somewhat less than the measurements (~10 cm"') at 
~700K for pyrene and coronene Poblin et al.|1995) . However, 
our "choice" to fit the observed, inherently asymmetric line pro- 
files of the out-of-plane bending modes with symmetric cal- 
culated profiles forces us to adopt a somewhat small intrinsic 
line profile. We will assess the effect of these assumptions on 
our fitting results later on. Lastly, we mention that the relative 
strength of the bands in the calculated spectrum will also de- 
pend on the internal energy (eg., temperature) of the emitting 
species and hence on the absorbed photon energy and the tem- 
perature cascade. However, over the limited wavelength range 
considered here, this effect is very small and we will here simply 
adopt a single absorbed UV photon energy (6.5 eV) characteris- 
ticfor the be nign conditions of the PDR in NGC 7023 Pobhn &| 
Mulas|2009) l. This corresponds to a peak temperature of ~900K 
for a 50 C-atom PAH. The database evaluates the temperature 
for each PAH and follows the temperature cascade consistently. 
Although our wavelength range is limited, we choose to use the 
temperature cascade to represented the most physically accurate 
approach. 



terms of abundance. Figure [8] shows a comparison between the 
PAH^ and PAH*^ contributions of the database fit (top), and the 
PAH^ and PAH'^ contributions in the BSS decomposition (bot- 
tom). The results of the database fits highlight various trends. 
First, except for the VSG spectrum, the overall contribution of 
large versus small PAHs seems very even. This may arise from 
the three times greater number of smaller PAHs in the database. 
Remarkably, although they represent only 25% of the database, 
the main features of each spectrum are generated by the larger 
PAHs while the continuum and less dominant features are re- 
produced primarily by the smaller PAHs. Second, both the fit 
of the observed raw spectra shown in Figure [8] and the vari- 
ous PAH charge state contributions shown in the right side of 
Figure [TJclearly shows that the dominant contributors to the 11.0 
fim band are PAH^ while the bulk of the 11.2 fim feature is due 
to PAH° species. This behavior, based on the spectra of hun- 
dreds of PAHs, confirms the early suggestion, based on a hand- 
ful of experimental PAH spectra, that the emission between 10.8 
and 11.1 fim can be used as a tracer of PAH cations (Hudgins 
|& Allamandola|1999[ ). It should also be noted that the BSS ex- 
tracted VSG signal has two strong absorption features at about 
1 1.0 and 12.6 jum. These probably arise from cross mixing as de- 
scribed earlier and are likely artificial (Appendix [A|. Therefore, 
the fit was also performed with a linear interpolation over these 
points. This fit produced almost identical results. 

Common concerns regarding the PAH Database fitting meth- 
ods are uniqueness and degeneracy. In the case of PAHs with a 
limited wavelength range, these concerns must be approached 
in a different way. This is further discussed and investigated in 
Appendix [B] These studies show that if the database is tasked to 
fit the BSS extracted PAH^ spectra with only PAH" species, this 
is not possible, further strengthening the attribution of the 11.0 
fim feature to PAH^. On the other hand, if the database attempts 
to fit the BSS extracted PAH" signal with only PAH+, a suitable 
fit is provided. However, without any limitations of the database 
fit, the reduced norm fit chooses predominately PAH^ to fit the 
BSS extracted PAH+ signal and PAH" to fit the BSS extracted 
pah" signal (Table [T}. This only shows that PAH cations have 
emission features that peak through the 10 - 15 fim range and it 
is important to employ more than one method to apply as many 
astronomical constraints as possible. 

6. Origins of AIB Variations in the 10-15 jum Range 

Investigating the spectra and examining the spatial distribution 
of the extracted emission components, we recognize four aspects 
that vary with spatial position: the [ 1 1 .0]/[ 1 1 .2] yum ratio, the red 
wing of the 11.2 jum feature, the precise peak position of the 11.2 
//m band, and the weak features at 12.0, 12.7, and 13.5 jim. In 
the following section we discuss these variations in the context 
of Signal 1, 2, and 3, which will now be referred to as PAH", 
PAH^, and VSGs, as emission feature carriers. 



5.2. Fit Results 

The results for the fits of the three extracted BSS component sig- 
nals are presented in Figure|7] broken down in categories of PAH 
size and charge. The small ( <50 C atoms) and the large (>50 C 
atoms) PAHs, as well as the cation, neutral, and anion species 
were separated out in order to judge each subgroups contribu- 
tion to the fit. A fit to the raw observed spectra at a typical po- 
sition was also made (Figure [8]l. The results are summarized in 
Table[T] It should be emphasized that the contribution of individ- 
ual species was quantified in terms of emission intensity, not in 



6.1. Ratio of 11.0 to 11.2 Features 

In most observations of the 10 - 15 //m range, the dominant 
11.2 //m feature appears with a dwarfed satellite feature at 
11.0 yum. This has been attribute d to the solo C-H out-of-plane 
bending modes of PAH+ (|Hudgins & Allamandola||1999[ 



van 



Diedenh oven et al.||2004[ |Hony et al.||2001[ [Bauschlicher et"aT 



i2008n200 9). We observe the highest [11.0]/[1 1.2] yum ratio clos- 
est to the source star, which is also where the abundance of PAH^ 
is greatest. A comparison of observed spectra from a PAH^ 
dominated region and a VSG dominated region are shown in 
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Fig. 7. The fit from NASA Ames PAH IR Spectroscopic Database (Bauschlicher et al. ] |2010| l of the three components extracted 
from BSS, PAH+ (top), PAH° (middle), and VSG (bottom). In the left column, we compare large PAH contiibution to small PAH 
contribution. In the right column we compare cation, neutral, and anion contribution to the fit. 
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Fig. 9. The original spectra from pixel positions (13, 15), where 
PAH^ are highly concentrated and (20, 10), where VSGs are 
abundant highlighting the extremes of the observed spectral vari- 
ations. 



Figure |9] It is important to recall that regardless of position in 
the PDR, the 11.2 fim PAH" feature significantly dominates the 
11. fim PAH^ feature. The separate contributions of PAH^ and 
PAH*^ to the 11.0 and 11.2 fim features respectively, agrees with 



previous spectroscopic laboratory and quantum chemical cal- 
culation studies by e.g. Hudgins & Allamandola ( 1999|l; Hony 
et al.'(2001');'van Diedenho ven et al.| ( |2004| ); |Bauschlicher et al. 
(2008) and Cami (2010). In our decomposition analysis, the 11.0 
jjm band is clearly associated with the PAH^ component (Signal 
2) based upon the strong 6.2 and 7.7 jim bands, while the 11.2 
fim band is attributed to the neutral component ( [Berne et al. 
|2007 ). Hence we suggest that the variation of the [11.0]/[11.2] 
fim ratio is due to a changing abundance of PAH+ to PAH". 



6.2. The 11.2 fim Red Wing 

As mentioned in the introduction, the 1 1.2 jjm band has an ob- 
served asymmetry with a varying red wing (Roche & Aitken 
1985a I. This wing has been attributed to anharmonicity or to dif- 
ferent species of PAHs with a shifted solo mode peak emission 
( |Pech et al.|2002[[van Diedenhoven et al.|2004f Observations of 
the 1 1 .2 fim feature show that the shape and peak position can 
vary. In their analysis of the skewed variations in the 11.2 fim 
profile, |van Diedenhoven et al. (j2004) empirically divided ob- 
servations of the feature into two categories: one group is char- 
acterized by a peak between 11.20 and 11.24 fim and a more 



7 



Rosenberg et al.: Variations of 10-15 yum AIBs of NGC 7023 



c 

Z) 

L_ 

o 

-t-l 
JD 

O 

X 



0) 
_N 

"5 



Raw Observations 



PAH Database Fit 
PAH^ 



PAH' 




BSS Method 




10 11 12 13 14 
wavelength[/i,m] 



11 12 13 14 
wavelength[/2m] 



11 12 13 14 
wavelength[yu,m] 



1 5 



Fig. 8. Comparison of the fitting results of the two methods, the NASA Ames PAH IR Spectroscopic Database (Bauschlicher et al. 
2010 1 (top) and BSS (bottom). The left panel depicts an observed spectra which is fit with both the database and the BSS extracted 
spectra. The middle and right panels are the cation (middle) and neutral (right) components obtained from the respective methods. 
The signals are normalized to the 1 1 .2 /urn peak. 



skewed red wing, while the other peaks at 1 1 .25 yum and is much 
more symmetric. In agreement with this, it has also been shown 
that PAH anions bands fall on the red side of the 11. 2 jum peak 
and could contribute to the red wing (Bauschli cher et al.|2008| l. 

By separating the observed spectra into component signals, 
we found that the main carrier of the 11.2 fim emission feature 
is pah". However, the observed spatial variations in the profile 
result in a PAH" source signal that is mainly symmetric with only 
a weak anharmonic red wing. The BSS analysis shows that the 
red wing is mainly due to a changing contribution of the VSG 
to the observed spectra; e.g. as the VSG signal becomes more 
prominent towards the outer edge of the PDR, its contribution 
to the observed spectra also increases (c.f. Figure |6|l. We note 
that in NGC 7023, position which are near to the exciting star 
are characterized by a more symmetric 11.2 fim profile while 
positions further away have a more pronounced red wing (cf.. 
Figure [9]l. 

The feature at 11.2 jim has been observed to shift peak posi- 
tion between 11.2 to 11.3 fim. Studying Figure |6] there are clear 
emission contributions from each component species throughout 
the 1 1 .0 to 1 1 .3 jum range. We propose that the variation of abun- 
dances of VSGs and PAHs in the observed spectra causes the 
shifting peak position of the 1 1.2 fim band. If there is a stronger 
contribution of PAHs in a certain region, the peak is observed to 
be blue-shifted. If the VSGs become more abundant, the peak is 



redshifted. This is in disagreement with [van Diedenhoven et aL| 
(2004) , where they notice a redshifted peak with a more sym- 
metric profile. 



6.3. The 12.0, 12.7, and 13.5 fim Features 



|Hony etaL] ( |2001| l assigned each emission feature to a differ- 
ent geometry and composition of PAH, depending on how many 
adjacent C-H groups are attached to the ring e.g., solo, du o, trio, 
and quartet modes of PAH" and PAH+ . The results by Hony et al. 
(2001) were further expanded to include compact and irregular 

which 



20091 



shaped large PAHs by Bauschlicher et al^( 2008 
are more astronomically relevant. [Bauschlicher et al.| ( |2009[ ) at- 
tributed the 11.3 - 12.3 fim band to the "duol" CH mode while 
the 12.5 - 13.2 fim region is attributed to "duo2" CHqop bands, 
the split of duo modes being caused by coupling to other bend- 
ing modes. The 13.5 fim feature has been attributed to the CHqop 
quartet mode of large irregular PAHs ( Bauschlicher et al.|2009 



Hon y et ar]|2001| l and can be used to place constraints on the 
edge structures of the emitting PAHs. Here we will place further 
astronomical constraints on the results of Hony et al. (2001 1 and 
IBauschhcher et al.| ( |2009l ). 
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6.3.1. The 12.0yum Feature 

The peak position of the "12.0 fim" feature varies from 1 1.8 yum 
to 12.0 fim (Figure |9]l. The 12.0 yum feature has been attributed 
to both the PAH° and PAH+ duo modes (|Hony et ali]|200r 



Through the separation of source signals, we have isolated the 
main 12.0 fj.m feature to the PAH*^. There is a feature that shares 
the profile of the 12.0 fim feature in the PAH^ spectrum but it is 
blue-shifted, peaking around 1 1.8 jim. 



6.3.2. The 12.7 yum Feature 

The 12.7 yum feature has been predominately attributed to the 
overlap of PAH^ duo and trio modes, but could not be def- 
initely attributed to eithe r PAH " or PAH+ ( |Hony et aL||2001 



Bauschlicher et al. 2008 2009 1. Examining the results of the 
signal separation in Figure |6| and Figure [8] two unique features 
at 12.7 yum are revealed, that of PAH+ and PAH". Although they 
share roughly the same peak position, the PAH" 12.7 /im feature 
is shifted to the red and is seen from 12.5 to 13.0 fim, while the 
PAH^ feature is blue shifted and asymmetric located between 
12.3 and 12.8 yum. We can attribute the variable blue wing of the 
12.7 fj.m feature to the changing abundance of PAH^. A peak 
shift and prominent asymmetry is seen in Figure |9] in Position 

I, located in the most concentrated area of PAH+. This feature 
is seen along with the increased 11.0 yum feature, a blue-shifted 
peak position of the 11.2 yum feature, a decreased red wing of the 

I I. 2 yum feature, and a "12.0 /im" feature peaking at 11.8 yum. 



6.3.3. The 13.5 yum Feature 

The observed spectra show a distinct 13.5 yum feature. In a study 
of Ml 7 it was suggested that this feature is coupled to the warm 
dust continuum ( Verstraete et al.|1996 l. Hony et al. (2001 1 fur- 
ther investigated this possibility and instead, attributed the 13.5 
yum feature to a quartet out-of-plane bending mode of PAH^ and 
PAH". Using BSS, we isolated this feature to PAH + (Signal 2) 
and pah" (Signal 1), in agreement with the results of Hony et al. 
( |2001| l, and likely decoupled from the warm dust continuum. 



6.4. Systematic Blue Shift witli Ionization 

We have attributed the 1 1.0 yum feature to PAH^, while the 1 1.2 
fim feature is attributed to PAH". In addition, the 11.8 yum fea- 
ture is attributed to PAH^, while the 12.0 yum feature is at- 
tributed to pah". We also identify the broad 12.7 yum band in 
both PAH^ and PAH", yet it appears to be a blend of features. 
The PAH-*- 12.7 yum band is also bluer than the PAH" 12.7 fim 
band. Specifically, the PAH^ broad 12.7 yum feature spans 12.3 
to 12.8 yum while the PAH" band stretches from 12.5 to 13.0 yum. 
Comparing the PAH^ and PAH" spectra, there is a systematic 
0.2 yum blue shift between the emitting bands. We do not observe 
this shift in the 13.5 yum band. PAH band shifts can occur due to 
temperature change in the emitting region, yet according to the 
model proposed by Pech et al. ( 2002 ), a 0.2 yum shift of the 11.2 
yum feature corresponds to a 650 K PAH temperature change. 
This PAH temperature change is too great to be observed within 
NGC 7023 NW, therefore it is unlikely that this band shift is due 
to a temperature change. Instead, we conclude that this shift is 
due to ionization, which modifies intrinsic emission properties of 
PAHs. Investigation on the exact origin of this shift is, however, 
beyond the scope of our paper. 



6.5. Other Possible Effects on the Shape of AIBs in the 10- 
15 fim Range 

Other effects and chemical properties have been reported to al- 
ter the shape of AIBs in the 10 - 15 yum range. Anharmonicity 



effects, as shown by e.g. Pech et al. 



( 2002| ) can modify the posi- 



tion and the symmetry of the 6.2 and 11.2 yum band and create the 
extended red wing in our observations. By means of DFT calcu- 
lations, [SiPAH]^ TT-complexes were also proposed by Joalland] 
et al. ( 2009 j ) to produce a splitting of the initial 11.2 yum PAH 
band into two bands at 1 1 .0 yum and 1 1 .4 yum due to the Si adsorp- 
tion on the PAH edge creates and a blue-shifted 6.2 fim band. 



We argue here (see Section 6.2 1, that the asymmetry of the 
11.2 fim feature is predominately due to the contribution of 
VSGs. Anharmonic effects are however still observed: the PAH" 
signal is not fully symmetric and displays a slight red wing, sug- 
gesting that anharmonicity effects are still important, but recall 
that most of the red wing is due to the varying abundance of the 
VSG component. Since [SiPAH]^ are expected to have a blue- 
shifted 6.2 yum band, we inspected the SL data but found no such 
signature. The splitting of the 11.2 feature is seen in Signal 2 
(PAH^), which is most concentrated in the regions near the star. 
Since the binding energy of [SiPAH]^ is about 2 eV, they should 
be destroyed easily the highly irradiated environment near the 
star Altogether, this suggests that compact PAH+ are a more 
natural explanation for the 11.0 fim feature, than [SiPAH]^ tt- 
complexes. 



6.6. Using the 11.0 and 1 1.2 fim features as tracers of 
ionization 

With the attribution of the 11.0 yum feature to PAH^ and the 
11.2 yum to pah", we can investigate the possibility of using 
this ratio to probe the ionization fraction of PAHs in the PDR. 
One of the classic methods to trace the PAH ionization fraction 



is the [6.2]/[11.3] fim integrated intensity ratio (e.g. Galliano 
|et al.|p008] l). There are other tracers of ionization such as the 
[7.7]/[11.3] fim and [8.6]/[11.3] fim ratios, but the 6.2, 7.7, and 
8.6 fim features include blended PAH^ and PAH" bands. As we 
show here, the 11.0 yum band is a purely cationic band and the 
11.2 yum band is purely neutral, increasing the accuracy of ion- 
ization fraction measurements. To demonstrate the reliability of 
the [1 1.0]/[1 1.2] yum ratio as an ionization indicator, we com- 
pare the [6.2]/[11.2] fim ratio to the [1 1.0]/[1 1.2] yum ratio us- 
ing the IRS-SL and IRS-SH observations of the NGC 7023-NW 
(Figure [T0| . In order to have an accurate measurement of the 
1 1.2 yum feature, without contamination from the 1 1.0 fim satel- 
lite feature, we compare the integrated intensity of the 6.2 yum 
feature from IRS-SL observations to the intensity of the 11.2 
fim feature using the high-resolution observations, since the 11.0 
was not resolved and separated in the IRS-SL observations. The 
maps were re-gridded using Montage so that each point of the 
SH map corresponds to the same spatial position on the SL map. 
Only the highest signal to noise data were used in this plot. 
For the 6.2 yum low-resolution map, we set a band integrated 
intensity threshold of 10 Wm^^sr ^ For the 11.0 yum high- 
resolution map we set a threshold of 10"^ Wm^^sr"' and the 1 1 .2 
fim high-resolution map has a threshold of 10"^ Wm^^sr The 
[6.2]/[11.2] vs [11.0]/[11.2] fim ratio in NGC 7023 is presented 



in Figure 10 The data reveal a clear correlation, validating the 



use of the [1 1.0]/[1 1.2] fim ratio as a PAH ionization indicator. 
The outliers in the upper left corner correlate to spectra where 
the thermal continuum from the source star is contaminating the 
linear continuum subtraction. For this reason, these points were 
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Fig. 10. The [6.2]/[11.2] fim ratio vs the [11.0]/[11.2] fim ratio in 
NGC 7023-NW. The 11.2 yum and 11.0 jim measurements were 
made using the IRS-SH observations while the 6.2 fim measure- 
ments were made using the IRS-SL observations. The circled 
data were not included in the fit (see text for details). The instru- 
mental error is comparable to the symbol size. 

not included in the linear fit. The linear fit has a high correla- 
tion coefficient of 0.95, from which an empirical relation can be 
derived: 



[U.Ofim] 
[ll.l^m] 



= 0.016 X 



[6.2fim] \ 
[ll.2fim]j 



(1) 



realize that there is a hidden issue here: the spectral differences 
in the 11-15 fim range - the broad and indistinct band in the VSG 
component versus the very distinct 11.2 and 12.7 modes of the 
PAHs - implies a more complicated evolutionary relationship 
between the VSGs and the PAHs than simple evaporation. 



8. Conclusion 

Applying a BSS method to observations from Spitzer's Infrared 
Spectrograph, Short-wavelength High-resolution mode, we un- 
covered 3 component signals in the PDR NGC 7023-NW. We 
found that each signal is most abundant in different regions of 
the PDR, depending on the radiation environment. We identified 
the three component signals as PAH cations, neutral PAHs, and 
VSGs. As the observed spectra suggest, the neutral PAHs domi- 
nate every region of the PDR, but are most heavily concentrated 
in between the PAH cations and VSGs. Both the spectra and 
spatial maps of each signal show high correlation to the results 
using Spitzers IRS-SL mode ( Berne et al.|2007 Berne and et al. 
2010), allowing us to use these results to verify our conclusions. 

To further explore the origin of the three resolved signals, 
we employed the NASA Ames PAH IR Spectroscopic Database. 
The fit shows that the component spectra resolved by BSS could 
be recreated by an appropriate combination of specific classes 
of PAH spectra from the database. Then, we used a database to 
fit an observed spectrum and grouped the individual molecules 
into charge class, then compared the spectra of the combined 
charge classes to the BSS extracted PAH^ and PAH" signals. The 
components were found to be very similar to the BSS extracted 
PAH+ and PAH°. 

Specific spectral properties are found for each population: 



7. Nature of the VSGs 

The BSS analysis identifies an independent broad component 
underneath the well-known 1 1 .2 and 12.7 yum bands. Earlier BSS 
studies over a wider wavelength range and the spatial distribu- 
tion of Signal 3 (Figure |5]l assigns this component to emission 
by VSGs, proposed to be PAH clusters. Early observations of 
the 1 1 .2 yum feature and its underlying emission support the sug- 
gestion that this broad underlying pedestal arises from a sepa- 
rate component (Roche & Aitken 1985b). The PAH Database 
analysis provides some further insight in the character of the 
carrier of this broad component. In this analysis, the 11-15 fim 
pedestal emission is due to a large number of individual com- 
ponents originating in a wide variety of molecular edge struc- 
tures (solo's, duo's, and trio's), which together blend in an in- 
distinct broad emission bump from 11-15 /urn. For this blend, 
the analysis selects relatively small species from the database. 
However, that is a selection effect. Small PAHs have, by neces- 
sity, a preponderance of corner structures. In contrast, calcula- 
tions for large PAHs have focused (for obvious reasons) on reg- 
ular structures with long straight edges and consequently strong 
1 1.2 yum bands and weak bands at longer wavelengths. We sur- 
mise that large irregular PAHs would equally fit the bill. The 
VSG component has been assigned to clusters of PAHs based 
upon an interpretation of the observed spatial distribution and the 
physical properties of clusters ( |Berne et al. 2007 RapacioU et al. 
2005||2006l). However, the spectroscopic properties of PAH clus- 



ters are presently unknown. While in general their spectra might 
be expected to resemble those of the constituent PAH molecules 
making up the cluster, we surmise that steric hindrance may af- 
fect the frequencies of the out-of -plane CH bending modes. We 



- We have attributed the 11.0 and 11.2 yum bands to cations and 
neutral species respectively. 

- We conclude that the variation of the [11.2]/[11.0] yum ratio 
depends on the relative abundances of PAH cations to neutral 
PAHs. 

- The extended red wing seen on the 11.2 yum feature is at- 
tributed to the increasing abundance of VSGs and the broad 
11.3 yum feature that is characteristic of this component. 

- The changing peak position of the 1 1.2 fim feature can also 
be explained by varying contributions from PAHs (blue shift) 
and VSGs (red shift). 

- The 12.0 yum feature is attributed to neutral PAHs while the 
11.8 yum feature is attributed to PAH cations, therefore, as the 
ionization mixture changes, the peak of this feature will shift 
accordingly. 

- Since the 13.5 yum feature is present in both PAH cations and 
neutral PAHs, but not existent in the VSG signal, where we 
see the continuum, we agree with Hony et al. (2001 1 that the 
13.5 yum signal is decoupled from the 15 fim continuum. 

By using the BSS method and the PAH Database fit, we 
arrived at the above conclusions. Each method has unique yet 
complementary strengths and weaknesses. The BSS method is 
blind, i.e. has no intrinsic assumptions about the emitting com- 
ponents, however since the statistical properties of the emit- 
ting components are unknown, the unmixing is not perfectly 
efficient. Additionally, the BSS method separates 3 mathemat- 
ically distinct signals, but gives no intuition about the molecular 
properties of these signals. The PAH Database allows for direct 
physical interpretation of the fit yet is biased towards smaller 
molecules and lacks spectral information for PAH clusters or 
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Other possible carriers of the VSG signal. Although both meth- 
ods suffer limitations, the strengths of one compensate for the 
weaknesses of the other Although the database fit may be de- 
generated in some cases, the interpretation of the database 
fit results, in terms of classes, is in agreement with the result of 
the BSS for PAH+ and PAH°. The VSG spectrum can only be 
obtained by BSS, and the database fit of this spectrum provide 
additional information on the possible chemical nature of this 
component. Both methods, i.e. BSS and Database Fitting, are 
powerful tools, but they must be used with an understanding of 
their limitations (described in details in Appendices A and B). 
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Appendix A: Exploring tlie Artifacts of BSS 

The efliciency of NMF is subject to two main limitations: 1) the 
possible non-unicity of solutions 2) the inaccuracy of the unmix- 
ing in the presence of noise. These two problems are the subject 
of intensive theoretical research in the field of signal processing 
(see e.g. Donoho & Stodden_2003 and Laurberg et al.^2008 ). If 
the statistical properties of the source matrix are known, then one 
can determine whether NMF will function properly. In a "real 
life" case like here however, since the matrix is not known a 
priori, we cannot verify these statistical properties. We therefore 
rely on some empirical considerations. Procuring the component 
spectra (Figure |3]l, the most prominent unmixing artifact is lo- 
cated at 11.0 jum and reveals itself with a sharp drop in Signal 
3. Some less discernible artifacts can be found at 11. 2 /im as an 
added peak of Signal 2 (seen in the grey envelope of Figure |3]l 
and at 12.7/im as a dip in intensity of Signal 3, roughly resem- 
bling an absorption feature. 

We suggest that the unmixing artifacts seen in some signals 
are compensated for by an increase or decrease of intensity, at 
the same wavelength, in other signals. This is most clearly seen 
at the 11.0 [ivs\ wavelength. The sharp drop to of Signal 3 is 
compensated by the weak satellite feature of Signal 1. Similarly, 
the absorption-type feature seen at 12.7 is compensated by a 
slightly increased intensity of Signal 2 at 12.7 /im. In order to 
better understand the unmixing efliciency, two tests were con- 
ceived: 



1. We have artificially mixed the database PAH*^ and PAH^ 
(Figure [8| signals with the NMF extracted VSG signal of 
[Berne et al. (2007) (Figure [6|. We mix these signals with a 
3 X 100 random matrix, but we make sure that the coeffi- 
cients correspond to realistic relative abundances (e.g. there 
are no regions with 100 % PAH^). We add synthetic white 
gaussian noise to the spectra to obtain a signal to noise ra- 
tio of about 20 (similar to our observations). In this way, we 
created 100 simulated spectra. We then apply NMF to this 



set of spectra. An example of the extracted signals is shown 
in Figure A. 1 There is a clear absorption band in the VSG 
spectrum, very similar to what is observed in the VSG spec- 
trum extracted from the observations (Figure |3]l. Lowering 
the noise does not improve these results. This implies that 
this artifact is caused by improper unmixing, and that this is 
an inherent problem of the NMF algorithm: the source spec- 
tra we are trying to extract most likely do not verify com- 
pletely the statistical properties that are needed for NMF to 
function efficiently ( jPonoho & Stodden|2003| l. We conclude 
that the absorption feature we observe in the VSG extracted 
spectrum (Figure [3]l, is caused by the improper unmixing of 

NMF 

In order to test why [Berne et al. ( 2007| did not see such an 
artifact in their VSG spectrum, the same BSS NMF algo- 
rithm was used to separate the IRS-SL data cube components 
but with only the 10-14 /vm wavelength range included. The 
result of this BSS is shown in Figure [ATS} where the 11.0 
/im absorption artifact is clearly visible. The fact that this 
artifact becomes present when only including the 10-14 /im 
spectral range implies a lack of information and constraints 
in this particular wavelength range. Rapacioli et al. ( 2005| l 
and ,Berne et al.| P007[ l show that the main features of the 
VSG spectrum are 1) a broad 7.7 //m band, 2) a continuum, 
dominant from 15 to 20 //m. Without these constraints in the 
reduced spectral range of IRS-SH, BSS has a difficulties ex- 
tracting enough information for this signal. We conclude that 
the artifacts especially in the VSG spectrum are mostly due 
to this reduced spectra range. 



Appendix B: Exploring the Limitations of the PAH 
Database Fitting 

To investigate the level of degeneracy of the database fit, we 
have performed several tests. Here, we present four fits with ex- 
plore this important issue. All fits follow the same values de- 
scribed in Section [sj a 6 cm ' FWHM, a 15 cm ' band red- 
shift, a Lorentzian profile, a 6 eV excitation energy with cascade 
emission, and molecules including only Carbon, Hydrogen, and 
Nitrogen. The fits are described below. 



1. 



The BSS extracted PAH^ spectra is fit with the database re- 
stricting it to only ionized species (top left of Figure ' 
The BSS extracted PAH" spectra is fit with the datab 



B.l 



stricting it to only neutral species (top right of Figure B.l i. 
The BSS extracted PAH^ spectra is fit with the database re- 



stricting it to only neutral species (bottom left of Figure B.l 
The BSS extracted PAH° spectra is fit with the database 
restricting it to only ionized species (bottom right of 
Figure B.l 



While the database fit of the BSS PAH+ signal with PAH+ 
cations is excel lent, when it is limited to neutral species (bot- 
tom left Figure B.l i, it is clear that the PAH^ spectrum cannot 
be recreated. However, the database fit of the BSS PAH*' species 
with cations does provide a suitable fit, except for one key re- 
gion: the 1 1.2 /im feature. The best fit of the database produces 
a two-pronged feature with neither prong peaking exactly at the 
11.2 /im band. However, it is possible that by varying the FWHM 
and band shift, which are somewhat arbitrary values, these bands 
could blend to recreate the 11.2 //m feature. This is due to the 
fact that various PAH cations peak throughout the mid-IR spec- 
trum. Specifically, PAH cations have emission features that peak 
throughout the 10.8 - 11.3 //m range. This only highlights again 
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the need for both approaches to apply as many astronomical 
constraints to the emitting species as possible, including wider 
spectral coverage when available. The BSS method allows us to 
spatially separate the three emitting components and therefore 
we create a physical picture in which Signal 2, located nearest to 
the star, represented the PAH^. In this case, the database supports 
this claim since we cannot fit this spectra with only neutral PAH 
species. We then move to Signal 1, which has a much different 
emission spectrum than the PAH^ and is located in the middle 
of the PDR with a much higher abundance than the PAH^. The 
database shows that this can be fit with either only neutral or only 
cation PAH species, but the physical intuition of BSS indicates 
that this is most likely a neutral PAH spectrum. 

Essentially, the database is a collection of Lorentzians that 
when restricted to a narrow band pass in which aromatic bands 
fall, and added together with varying intensity, can likely fit 
any type of aromatic-rich materials. Especially in the case of 
cations, the fit is degenerate (Figure [BTT] i. Yet when the fit is per- 
formed with the full range of the database, the highest solu- 
tion is strikingly similar to the BSS extracted solution (Figure[8]l. 
Although our tests show that the fits are poorly constrained for 
individual PAHs, this is not the case when considering PAH 
sub populations, which trace properties such as size and charge. 
Thus, the astronomical mid-lR bands studied here are assigned 
to PAH molecules as a class, not to individual PAH molecules. 
The literature is rich in examples where it is shown that the mid- 
IR PAH emission features can be categorized by charge, size. 



shape and composition (Peelers et al. 


2004 Bauschlicher et al. 


2008 2009 Boersmaetal.|2010l. W 


len a single species is re- 



moved from the pool of species that can be fit, another species 
from the same subgroup will simply replace it. Therefore, a mid- 
IR fit made with the database is expected not to provide infor- 
mation about the presence of an individual PAH molecule, but 



about the subclasses of PAHs involved and the response of that 
subclass to the local astronomical environment. 

Aside from charge, the database also probes the subclass of 
size. Although we have a limited wavelength range, the main 
features are best reproduced by larger PAH molecules. Similarly 
to the case of charge, we can treat PAH size as another sub pop- 
ulation with a set of characteristic properties. In Figure |B.2| we 
use the database to fit the observed spectra from Figure fej Here, 
besides size the only constraint applied is including only C, H, 
and N atoms. The parameters are equal to those described in 
Section 5.1. The top panel displays a fit of the observed spec- 
trum using only PAHs with 30 or more C atoms. The bottom 
panel is a fit of the observed spectrum using PAHs with 30 or 
less C atoms. There are 255 PAH species available to fit C < 30, 
whereas there are 237 PAH species available to fit C > 30, once 
the other constraints are taken into account. Even though there 
are less PAH species included in the fit of large PAHs, it is a su- 
perior fit to the small PAHs and has a lower reduced norm value. 
From this limited wavelength range, we can already deduce that 
large PAHs are responsible for the main bands, however to have 
a stronger constraint on size, an analysis using the full 3-15 jum 
range in necessary. 
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Fig. A.2. The result of BSS using the NMF method on the IRS-SL data from[ 

10-14 fim. BSS extracted PAHh- in blue, extracted PAH° in green, and extracted VSG in red. 
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Fig.B.l. PAH Database fit results of fitting the BSS extracted PAH^ (left) and PAH" (right) signals using only cations (blue) or only 
neutral species (green). The black line represents the original BSS extracted spectra while the fit is represented in color. 



14 



Rosenberg et al.: Variations of 10-15 yum AIBs of NGC 7023 




Fig. B.2. PAH Database fit results of fitting an observed spectrum limiting the species to large PAHs (C > 30) (top) and small PAHs 
(C < 30) (bottom). The black line represents the original observed spectrum, where the fit is represented in color. 
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